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Abstract. We have constructed the leading order hyperon-nucleon potential in a chiral Effective Field Theory ap- 
proach. The chiral potential consists of one-pseudoscalar-meson exchanges and non-derivative four-baryon contact 
terms. The hyperon-nucleon interactions are derived using SU(3) symmetry, the nucleon-nucleon interaction is not 
considered explicitly since it can not be described well with a leading order chiral Effective Field Theory. We solve a 
regularized Lippmann-Schwinger equation and show that a good description of the available hyperon-nucleon scatter- 
ing data is possible with five free parameters. The chiral potential can be used as further input for hypernucleus and 
hypernuclear matter calculations. Preliminary hypertriton calculations yielded the correct hypertriton binding energy. 
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1 Introduction 

Since the work of Weinberg [ 1 , 2 1 , the derivation of nuclear in- 
teractions from chiral Effective Field Theory (EFT) has been 
discussed extensively in the literature. For reviews we refer to 
00. The main advantages of this scheme are the possibilities 
to derive two- and three-nucleon forces as well as external cur- 
rent operators in a consistent way and to improve calculations 
systematically by going to higher orders in the power counting. 

Recently the nucleon-nucleon (NN) interaction has been 
described to a high precision using chiral EFT ||5][6). In |6], 
the power counting is applied to the NN potential, as origi- 
nally proposed in [1 2|. The NN potential consists of pion- 
exchanges and a series of contact interactions with an increas- 
ing number of derivatives to parameterize the shorter ranged 
part of the NN force. A regularized Lippmann-Schwinger equa- 
tion is solved to calculate observable quantities. Note that in 
contrast to the original Weinberg scheme, the effective poten- 
tial is made explicitely energy-independent as it is important 
for applications in few-nucleon systems (for details, see Q). 

The hyperon-nucleon (YN) interaction has not been inves- 
tigated using EFT as extensively as the NN interaction. Hy- 
peron and nucleon mass shifts in nuclear matter, using chiral 
perturbation theory, have been studied in [8|. These authors 
used a chiral interaction containing four-baryon contact terms 
and pseudoscalar-meson exchanges. Recently, the hypertriton 
and Ad scattering were investigated in the framework of an 
EFT with contact interactions [9|. Some aspects of strong AN 
scattering in effective field theory and its relation to various 
formulations of lattice QCD are discussed in iflOl . Korpa et 
al. ifTTl performed a next-to-leading order (NLO) EFT analysis 
of YN scattering and hyperon mass shifts in nuclear matter. 
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Their tree-level amplitude contains four-baryon contact terms; 
pseudoscalar-meson exchanges were not considered explicitly, 
but SU(3) breaking by meson masses was modeled by incorpo- 
rating dimension two terms coming from one-pion exchange. 
The full scattering amplitude was calculated using the Kaplan- 
Savage-Wise resummation scheme fl2l . The YN scattering 
data were described successfully for laboratory momenta be- 
low 200 MeV, using 12 free parameters. 

In this contribution we show the results for the recently con- 
structed chiral EFT for the YN system fPJl . In this work we 
apply the scheme used in [6 1 to the YN interaction. Analogous 
to the NN potential, at leading order (LO) in the power count- 
ing, the YN potential consists of pseudoscalar-meson (Gold- 
stone boson) exchanges and four-baryon contact terms, related 
via SU(3) symmetry. We solve a regularized coupled channels 
Lippmann-Schwinger equation for the LO YN potential and fit 
to the low-energy YN cross sections, which are dominated by 
S'-waves. We remark that our approach is quite different from 

nm 



2 The effective potential 



In this section, we discuss in some detail the effective chi- 
ral YN potential at leading order in the (modified) Weinberg 
power counting. This power counting is briefly recalled first. 
Then, we present the minimal set of non-derivative four-baryon 
interactions and show the formulae for the one-Goldstone-boson- 
exchange contributions. 
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2.1 Power counting 

We apply the power counting to the effective hyperon-nucleon 
potential V e g which is then injected into a regularized Lippmann- 
Schwinger equation to generate the bound and scattering states. 
The various terms in the effective potential are ordered accord- 
ing to 



V eS = Ves(Q, g,(J.)=J2 & MQ/H 9) 



(1) 



where Q is the soft scale (either a baryon three-momentum, a 
Goldstone boson four-momentum or a Goldstone boson mass), 
g is a generic symbol for the pertinent low-energy constants, [i 
a regularization scale, V„ is a function of order one, and v > 
is the chiral power. It can be expressed as 



B + 2L 



Vi A: 



(2) 



with B the number of incoming (outgoing) baryon fields, L 
counts the number of Goldstone boson loops, and vi is the num- 
ber of vertices with dimension Ai. The vertex dimension is ex- 
pressed in terms of derivatives (or Goldstone boson masses) di 
and the number of internal baryon fields bi at the vertex un- 
der consideration. The LO potential is given by v = 0, with 
B = 2, L = and A.- L = 0. Using eq. (O it is easy to see 
that this condition is fulfilled for two types of interactions - 
a) non-derivative four-baryon contact terms with bi = 4 and 
di = and b) one-meson exchange diagrams with the leading 
meson-baryon derivative vertices allowed by chiral symmetry 
(bi = 2, di = 1). At LO, the effective potential is entirely given 
by these two types of contributions, which will be discussed in 
more detail in the following sections. 

2.2 The four-baryon contact terms 

The LO contact term for the A A interactions is given by fl]0 
C = C, {NTiN) {NTiN) , (3) 
where are the usual elements of the Clifford algebra |[T4| 

A = i , r 2 = y , r 3 = a» v , r 4 = 7 " 7s , r 5 = 75 . (4) 

Considering the large components of the nucleon spinors only, 
the LO contact term, eq. (0), becomes 

C = (yV) (vV) ~ -f (^W) {f^tp) i ( 5 ) 

where ip are the large components of the Dirac spinor and Cs 
and Ct are constants that need to be determined by fitting to the 
experimental data. In the case of YN interactions we will con- 
sider a similar but SU(3)-invariant coupling. The LO contact 
terms for the octet baryon-baryon interactions, that are Hermi- 
tian and invariant under Lorentz transformations, are given by 
the SU(3) invariants, 

C 1 = C\ (B a B b (r t B) b (r t B) a ) , 

C 2 =C?(B a (r i B) a B b (r i B) b ) , 

£ 3 = Cf (B a (r t B) a ) (B b {r t B) b ) . (6) 




•A A 




-N N 



Fig. 1. Lowest order contact terms for hyperon-nucleon interactions 



Here a and b denote the Dirac indices of the particles, B is the 
usual irreducible octet representation of SU(3) given by 



B = 




(7) 



and the brackets denote taking the trace in the three-dimensional 
flavor space. The Clifford algebra elements are here actually 
diagonal 3 x 3-matrices. The LO YN contact terms given by 
these interactions are shown diagrammatically in fig. [T] The 
partial wave potentials are, up to an overall factor of 4ir, 

Viso = 2 (Cf — 3Cf.) + 2 (C§ - 3C|.) , 
= 2(Ci + C!)+2(Cf + Cf), 

Visa — g (Cs _ 3C T ) + - (C| - 3Cf.) + 2 (Cj — 3Cf.) , 

v 3 A si = - (c£ + ci) + (C| + C 2 T ) + 2 (cf + C 3 T ) , 

v iso 



V 3S1 



2 (Cf - 3C|) + 2 (Cg - 3(7|.) , 
-2(<? 



I + C T ) 



2 (Cf + C 3 



T7 > 



^iso 7 — 2 ~ ^Ct) _ (Cs — 3Cy) + 2 (C| 3Cf.) , 
^3sf = | + C£) + (C| + Cf ) + 2 (Cf + C 3 ) , 
^iso = 2 (^s ~ 3Cy) — (Cf — 3C|.) , 



[cl + c 2 T ) 



(8) 



for the A A -» A A, ylA -> AA, ZW -> ZW (1=3/2), 
ZW -> £JV (1=1/2) and AN -> ZW channels respectively. 

The six contact terms, C\, C T , C|, C T , Cg, C|., are de- 
fined analogously to eq.|5]and enter the A A and YN potentials 
in only 5 different combinations. The sixth combination occurs 
only in the AA, EN and US channels. This is equivalent to 
the fact that only 5 of the {8} x {8} = {27} + {10} + {10*} + 
{8} s + {8} a + {1} representations are relevant for AA and 
YN interactions. So, 5 contact terms need to be determined by 
a fit to the experimental data. Since the AA data can not be 
described well with a LO EFT, see lfTl [T5l . we will not consider 
the AA interaction explicitly. Therefore, we consider the YN 
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partial wave potentials 

\rAA _ rii 
V 1S0 — °1 



.1.1 

SO' 



V 1S0 — °1S0> 



AA 
SO 



8Cl 



SE 
SO' 



\rAA _ f~iAA 
V 3S1 — °3S1' 



V 3S1 ~ ^3S1 ' 



\rSS _ r<AA 
V 3S1 — °3S1' 



(9) 



\rAE _ o (r<AA _ r< ss \ V AS — nAS 
V 1S0 ~ ° V^ISO ^lSO/ ' V 3S1 — °3S1 - 

We have chosen to search for C^, C^g X , Cf^, Cfgf , and 
in the fitting procedure. The other partial wave potentials 
are then fixed by SU(3) symmetry. 



2.3 One pseudoscalar-meson exchange 

Since they are discussed extensively in the literature, we show 
only briefly the results for one-pseudoscalar-meson exchanges. 
The spin-space part of the LO one-pseudoscalar-meson-exchang 
potential is, similar to the static one-pion-exchange potential in 
(recoil and relativistic corrections give higher order contri- 
butions), 



V 



-f. 



B 1 B'pfB 2 B'P 



(<r i • k) ( ffa • k) 

k 2 + nil, 



(10) 



where rap is the mass of the exchanged pseudoscalar meson. 
We defined the transferred and average momentum, k and q, 
in terms of the final and initial center-of-mass (cm.) momenta 
of the baryons, p' and p, as k = p' p and q = (p' + p)/2. 
The isospin part of the interaction and the SU(3)-invariant cou- 
plings Jbbp are given e.g. in Q6:|. The physical r\ was identi- 
fied with the octet r\ (r\%) and its physical mass was used. 



3 Scattering equation and observables 

The calculations are done in momentum space, the scattering 
equation we solve for the YN system is the (nonrelativistic) 
Lippmann-Schwinger equation. The Lippmann-Schwinger equa- 
tion is solved in the particle basis, in order to incorporate the 
correct physical thresholds and the Coulomb interaction in the 
charged channels. The Coulomb interaction is taken into ac- 
count according to the method originally introduced by Vin- 
cent and Phatak [17] (see also [18|). We have used relativistic 
kinematics for relating the laboratory energy 7] a b of the hyper- 
ons to the cm. momentum. Although we solve the Lippmann- 
Schwinger equation in the particle basis, the strong potential is 
calculated in the isospin basis. The potential in the Lippmann- 
Schwinger equation is cut off with the regulator function, 



(11) 



in order to remove high-energy components of the baryon and 
pseudoscalar meson fields. The total cross sections are found 
by simply integrating the differential cross sections, except for 
the E + p — » E + p and E~p — » channels. For those chan- 
nels the experimental total cross sections were obtained via 
El 

cos 6> max - cos min J cos emin d cos 
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Fig. 2. "Total" cross section a (as defined in eq. d!2b ) as a function 
of piab- The experimental cross sections in a are taken from refs. [20| 
(open squares) and (2JJ (filled circles), in b from refs. |22| (filled 
circles) and 1 23 1 (open squares) and in c,d from 1 19 1. The shaded band 
is the Jiilich chiral EFT for A = 550, 700 MeV, the dashed curve is 
the Jiilich '04 model [24], and the solid curve is the Nijmegen NSC97f 
model [ 161. 



for various values of cos # m i n and cos (9 max . Following [ 16 1, we 
usecos# m i n = — 0.5 and cos max = 0.5 in our calculations for 
the S + p — > S + p and E~p — > S~p cross sections, in order to 
stay as close as possible to the experimental procedure. 



4 Results and discussion 

For the fitting procedure we consider the empirical low-energy 
total cross sections shown in figs. |2^,c, and d and [3^ and b, 
and the inelastic capture ratio at rest (29), in total 35 YN data. 
These data have also been used in 02411161 . The higher energy 
total cross sections and differential cross sections are then pre- 
dictions of the LO chiral EFT, which contains five free param- 
eters. The fits are done for fixed values of the cut-off mass (yl) 
and of a, the pseudoscalar FJ (F + D) ratio. 

The five LECs in eq. (0, were varied during the parameter 
search to the set of 35 low-energy YN data. The values ob- 
tained in the fitting procedure for cut-off values between 550 
and 700 MeV, are listed in table Q] We remark that the range 
of cut-off values is similar to the range in the NN case, e.g. 
[30,31 1. The range is limited from below by the mass of the 
pseudoscalar mesons. Since we do a LO calculation we do not 
expect a large plateau (i.e. a practically stable x 2 f° r varying 
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Fig. 3. As in fig. [2] but now the experimental cross sections in a,b are 
taken from refs. 1251 and in c,d from 1261 . 

Table 1. The YN S'-wave contact terms for various cut-offs. The 
LECs are in 10 4 GeV -2 ; A is in MeV. \ 2 is tne total cn i squared 
for 35 YN data. 



.1 


550 


600 


650 


700 


n AA 
'-'ISO 
r<AA 
3S1 

n ss 

°1S0 

n ss 

°3S1 
n AS 


-.0466 
-.0222 
-.0766 
.2336 
-.0016 


-.0403 
-.0163 
-.0763 
.2391 
-.0019 


-.0322 
-.0097 
-.0757 
.2392 
.0000 


-.0304 
-.0022 
-.0744 
.2501 
.0035 


x 2 


29.6 


28.3 


30.3 


34.6 



A). Note that recently the cut-off dependence of EFT calcula- 
tions has been discussed extensively in the literature, see e.g. 

emu. 

The fits were first done for the cut-off mass A = 600 MeV. 
We remark that the AN .S'-wave scattering lengths resulting 
for that cut-off were then kept fixed in the subsequent fits for 
the other cut-off values. We did this because the AN scatter- 
ing lengths are not well determined by the scattering data. As 
a matter of fact, not even the relative magnitude of the AN 
triplet and singlet interaction can be constrained from the YN 
data, but their strengths play an important role for the hyper- 
triton binding energy [35|. Contrary to the NN case, see, e.g. 
1 15 1, the contact terms are in general not determined by a spe- 
cific phase shift, because of the coupled particle channels in the 
YN interaction. Furthermore, due to the limited accuracy and 



Fig. 4. Differential cross section da jd cos 9 as a function of cos 9, 
where 9 is the cm. scattering angle, at various values of pi a b (MeV/c). 
The experimental differential cross sections in a,b are taken from 1 19 1 
and in c,d from 1251 . Same description of curves as in fig. [2] 



incompleteness of the YN scattering data there are no unique 
partial wave analyses. Therefore we have fitted the chiral EFT 
directly to the cross sections. A good description of the consid- 
ered YN scattering data has been obtained in the considered 
cut-off region, as can be seen in table[T]and figs.|2^,c,d and[3^,b. 
In these figures the shaded band represents the results of the 
chiral EFT in the considered cut-off region. In this low-energy 
regime the cross sections are mainly given by the S'-wave con- 
tribution, except for for the AN — ► EN cross section where 
the 3 Di(AN) 3 Si(SN) transition provides the main con- 
tribution. Still all partial waves with total angular momentum 
J < 2 were included in the computation of the observables. 
The Ap cross section shows a clear cusp at the S + n threshold, 
peaking at 60 mb. Fig. [2J3 shows that the predicted Ap cross 
section at higher energies is too large, which is related to the 
problem that some LO phase shifts are too large at higher en- 
ergies. Note that this was also the case for the TV TV interaction 
021. In a NLO calculation this problem will probably vanish. 
The differential cross sections at low energies, which have not 
been taken into account in the fitting procedure, are predicted 
well, see fig.|4] The results of the chiral EFT are also in good 
agreement with the scattering data at higher energy, the older 
ones in figs. [3};,d as well as the more recent scattering data in 

fig.H] 

We have, so far, used the SU(6) value for the pseudoscalar 
F/{F + D) ratio; a — 0.4. We studied the dependence on 



H. Polinder: FZJ-IKP-TH-2006-42 Hyperon-nucleon interactions in effective field theory 



5 




-0.5 0.5 
cos 9 



0.5 0.5 
cos 9 



100 



J2 

E 




300 400 500 600 700 800 
p lab (MeV/c) 

Fig. 5. Recent YN data, a and b: differential cross section da jd cos 8 
as a function of cos 6, where 6 is the cm. scattering angle, at various 
values of p\^ (MeV/c). The experimental differential cross sections 
are from 1271 and (28 1, respectively, c: "total" cross section a as a 
function of piat>- The experimental cross sections are from 1 28 1. Same 
description of curves as in fig. [2] 



this parameter by varying it within a range of 10 percent. After 
refitting the contact terms we basically found an equally good 
description of the empirical data. Therefore, we keep a to its 
SU(6) value. An uncertainty in our calculation is the value of 
the r] coupling, since we identified the physical -q with the octet 
rj. Therefore, we varied the r\ coupling between zero and its 
octet value, but we found very little influence on the description 
of the data (in fact, inclusion of the r\ leads to a better plateau 
in the cut-off range considered). Also baryon mass differences 
squared in the propagator in eq. (flOl l are not consistent with 
our power counting and were not considered, but we studied 
their influence on the description of the data. We found that the 
quality of the description does not depend on these terms. 

The S'-wave phase shifts for Ap and U + p are shown in fig. 
[6] As mentioned before, the limited accuracy of the YN scat- 
tering data does not allow for a unique phase shift analysis. 
This explains why the chiral EFT phase shifts are quite differ- 
ent from the phase shifts of the models presented in refs. 11241 
[161 . The predictions of the latter models also differ between 
each other in many partial waves. In both the Ap and E + p 1 Sq 
partial waves, the LO chiral EFT phase shifts are much larger 
at higher energies than the phases from 112411161 . We emphasize 
that the empirical data, considered in the fitting procedure, are 
at lower energies. Also for the NN interaction in leading order 




o 
-a 
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Fig. 6. The Ap and E+p S -wave phase shifts 8 as a function of pi a b- 
Same description of curves as in fig. [2] 



these partial waves were much larger than the Nijmegen phase 
shift analysis, see 03). It is expected that this problem for the 
YN interaction can be solved by the derivative contact terms 
in a NLO calculation, just like in the NN case. Our 3 Si S + p 
phase shift is repulsive like in 1 24 ] , but contrary to lfl6ll . 

The Ap and S + p scattering lengths are listed in table [2] 
together with the corresponding hypertriton binding energies 
(preliminary results of YNN Faddeev calculations from 1 36 1). 
Our singlet S + p scattering length is about half as large as the 
values found in 112411 1 61 . Similar to those models and other YN 
interactions, the value of the triplet U + p scattering length is 
rather small. Contrary to [16], but similar to IT241 we found re- 
pulsion in this partial wave. The magnitudes of the Ap scat- 



Table 2. The predicted YN singlet and triplet scattering lengths (in 
fm) and hypertriton binding energy, Eb (in MeV). Experimentally the 
deuteron binding energy is —2.224 MeV and the hypertriton binding 
energy is —2.354(50) MeV. The binding energies for the hypertriton, 
1361 , are calculated using the Idaho-N3LO NN potential 0. 



A 


550 


600 


650 


700 


af p 
af P 


-1.90 
-1.22 


-1.91 
-1.23 


-1.91 
-1.23 


-1.91 
-1.23 


S+p 


-2.24 
0.70 


-2.32 
0.65 


-2.36 
0.60 


-2.29 
0.56 


E B 


-2.35 


-2.34 


-2.34 


-2.36 
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tering lengths are smaller than the corresponding values of the 
Nijmegen NSC97f and Jiilich '04 models EMED, which is also 
reflected in the small Ap cross section near threshold, see fig. 
12^. The mentioned models lead to a bound hypertriton 03611371 . 
Although our Ap scattering lengths differ significantly from 
those of [24, 16], the YN interaction based on chiral EFT also 
yields a correctly bound hypertriton, see table |2] Preliminary 
results for the four-body hypernuclei and ^He, see [38], 
show that the chiral EFT predicts reasonable A separation en- 
ergies for \H, but the charge dependence of the A separation 
energies is not reproduced. 

Finally, our findings show that the chiral EFT scheme, ap- 
plied in ref. J6) to the NN interaction, also works well for the 
YN interaction. In the future it will be interesting to perform 
a combined NN and YN study in chiral EFT, starting with 
a NLO calculation. Also an SU(3) extension to the hyperon- 
hyperon (YY) sector is of interest. Work in this direction is in 
progress. 

I thank Johann Haidenbaur and Ulf-G. MeiBner for collaborating on 
this work. Also I am very grateful to Andreas Nogga for providing me 
with the hypernuclei results. 
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